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The Dyson hierarchical one-dimensional Ising model of parameter a > contains long-ranged 
ferromagnetic couplings decaying as l/r 1+CT in terms of the distance r. We study the stochastic 
dynamics near zero-temperature via the Real Space Renormalization introduced in our previous 
work (C. Monthus and T. Garel, arXiv: 1212.0643 ) in order to compute explicitly the equilibrium 
time teq(L) as a function of the system size L. For a < 1 where the static critical temperature for the 
ferromagnetic transition is finite T c > 0, we obtain that dynamical barriers grow as the power-law: 

lnt eq (L) ~ /3 ( 3(2^-° -i) ) ^ 1-<T - F° r f = 1 where the static critical temperature vanishes T c = 0, 

we obtain that dynamical barriers grow logarithmically as : lnt eq (L) ~ [f] ( gf 9 2 ) ~ l] hiL. We also 
compute finite contributions to the dynamical barriers that can depend on the choice of transition 
rates satisfying detailed balance. 

I. INTRODUCTION 

The Dyson hierarchical Ising model [l| has been introduced as a ferromagnetic model where the partition function 
could be analyzed via exact renormalization. The hierarchical couplings are chosen to mimic effective long-range 
power-law couplings J(r) ~ l/r 1+IJ in one dimension so that a phase transition with finite T c is possible in the region 
< a < 1. This type of hierarchical model has thus attracted a great interest in statistical physics, both among 
mathematicians [2H5j and among physicists Note that Dyson hierarchical models have been also introduced in 

the field of quenched disordered models, in particular for random fields [13, El) spin-glasses (12-14]. and for Anderson 
localization |15l422l |. 

In the present paper, we consider the Dyson hierarchical ferroma gnet ic Ising model and we focus on the stochastic 
dynamics satisfying detailed-balance, such as the Glauber dynamics [23j . We do not consider the coarsening dynamics 
starting from a random initial condition (see the review [24j |) but focus instead on the equilibrium time t eq near zero 
temperature, i.e. the time needed to go from one ground state (where all spins take the value +1) to the opposite 
ground state (where all spins take the value —1). In a previous work, we have introduce a real-space renormalization 
procedure to determine this equilibrium time t eq [25| as a function of the system size. Here we solve the corresponding 
RG flow for the Dyson hierarchical Ising model. 

The paper is organized as follows. In section UH we recall the important properties of the equilibrium of the 
Dyson hierarchical Ising model. In section IIII1 we study the stochastic dynamics satisfying detailed-balance via the 
real space renormalization procedure introduced in |25j [. The explicit solutions of the renormalization flow for the 
'simple' dynamics and for the Glauber dynamics are presented in sections IIVI and IVl respectively. Our conclusions are 
summarized in section fVTl Finally in Appendix [XI we discuss the link between dynamical barriers and the energy-cost 
of a single domain wall. 

II. REMINDER ON THE STATICS OF THE DYSON HIERARCHICAL ISING MODEL 

A. Definition of the Dyson hierarchical Ising model 

The Dyson hierarchical Ising model is a model of 2 N classical spins Si = ±1 where each configuration C = 
(Si, S 2 , ..S 2 n) nas f° r energy 
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where the positive couplings J n depend on n via 

J n = J 2^ n (2) 

To make the link with the physics of long-range one-dimensional models, it is convenient to consider that the sites 
i of the Dyson model are displayed on a one-dimensional lattice, with a lattice spacing unity. Then the site i = 1 is 
coupled via the coupling J n /(2 n ) 2 to each spin of index 2"~ 1 < i < 2™. At the scaling level, the hierarchical Dyson 
model is thus somewhat equivalent to the following power-law dependence in the real-space distance L n = 2™ 



rff(T ) - J " - J ° - J ° (1) 
1 n> ~ (2 n ) 2 ~~ 2™( 1 + CT ) ~~ L 1+a 



The parameter a is thus the important parameter of the model. 
The ground state energy where all spins have the same sign 

U^n — U 2 n(Si — 1, 5*2N =1) = — Jq2 1 — J\2 2 — J 2 2^ ^ — ... — Jn—i 

= -J 2 N - 1 — — = -^L N — ^- (4) 
1 - 2- a 2 1 - 2- a y ' 

is extensive in the number Ljv = 2 N of spins in the region 

a > (5) 

The energy cost of the configuration where the first 2 1 spins arc (—1), whereas the other 2 1 spins arc (+1) 

Ug" 1 ^ - U? N S = 2J N -t = 2J 2^ N ^ = J 2°L^° (6) 

grows with the distance Ln — 2 N for a < 1, remains constant for a = 1, and decays for a > 1. It has been shown [l[ 
that the critical temperature T c for the ferromagnetic transition is finite for a < 1 only 

T c (a < 1) > 

T c (a > 1) = (7) 

For a > 1 , note that Eq [5] is not very physical since the energy cost of a single domain wall in the middle of the 
sample becomes exponentially small instead of remaining finite in the usual one-dimensional ferromagnetic chain. So, 
in the following, we will consider only the region 

< a < 1 (8) 



B. Real-space renormalization for the equilibrium near zero-temperature 



In this article, we will focus on the regime 'near zero temperature' 

T< Jo 



(9) 



where the partition function is dominated by the two ferromagnetic ground states. Then the real-space renormalization 
of the partition function becomes very simple. The consecutive spins (S21-I) S 2 i) that are coupled by Jq are grouped 



into a single renormalized spin = 



S 2i -i+S 2 i 



corresponding to the ferromagnetic cluster of the two spins 



|£f = 1 > = I-Sm-i = 1 > \S 2i = 1 > 
IS? = -1 > = \S2i-1 = -1 > \S 2i = -1 > 



(10) 



In terms of these 2 N 1 renormalized spins for £ = 1,2, .., 2^ 1 , the energy of Eq. [T]becomes 

U 2 N (Si, S 2 n) = — Jq2 1 - Ji \S1S2 + S3S4 +...]- J; 
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J. 



N-l 



OR 



2 N- 2 



S^N -2 H~ ■•• H~ — 1 



2 N- 2 



Jo2 



N-l 



+ 2 1 a U 2 N-i (Si , S^n-i) 



(11) 



3 



So if one considers the partition function where the ratio U/T enters, one obtains that Eq. [TT] corresponds to the 
following renormalization of the temperature upon the elimination of the lowest generation 

T R = 2°- 1 T (12) 

This is in agreement with Eq. [7] : for a < 1, Tr flows towards the attractive fixed point T = 0; for a > 1, Tr flows 
away from the unstable fixed point T c = 0. 

To prepare the following sections, it is also convenient to describe the renormalization of the local field. In the 
initial model of Eq. Q] the local field on spin Si is defined as 

Bi = _ du 2 HSu s ..,s 2N ) = ^ JijSj (13) 

i.e. for instance for the two first spins i = 1, 2 



B\ — J0S2 + "7T 



■h 1 S5 + Sq + S7 + Sg \ t Jn-i ( S 2 n-x + i + ... + S 2 ^ 



2 V 2 J 2 2 V 2 2 / 2^-! V 2^-! 

Ji / 1S3 + Si\ J 2 f S 5 + S 6 + S 7 + S$\ Jn-i ( 5* 2 n-i +1 + ... + 5 2 jv \ 
S 2 = J 5 1 + y ^^— - j +... + __ ^ __ j (14) 

After the renormalization where (Si,^) have been grouped into the renormalized spin = (Si + S2V2 fEq [TO]) , 
the renormalized local field B^ on 5*^ reads 



B* =Ja 



03 + S4 \ f + D6 + S7 + Ss \ Jn-1 ( SoN-i+i + ■•• + S. 



2 l - a J 



2 J 2 \ 2 2 / 2 W - 2 V 2^-! 



(15) 



i.e. it is renormalized by the same factor 2 a as the couplings as it should. 

III. REAL SPACE RENORMALIZATION FOR THE DYNAMICS OF THE DYSON MODEL 

A. Stochastic single-spin-flip dynamics satisfying detailed-balance 

We consider the stochastic dynamics generated by the Master Equation for the probability -Pt(C) to be in configu- 
ration C at time t 

= J2 p t (O w (C -+c)-p t (c) Y, w ( c c ') ( 16 ) 

c c 
where the transition rates W (C — > C) satisfy the detailed balance condition 

e -W{c) w (c _> c /j = e -mc') W (c' C) (17) 
For the Ising model of Eq [l] where the energy is of the form 

U{C) = -J2JiiS i S j (18) 

i<j 

with the local fields (Eg [13]) 

it is natural to consider the following single-spin-flip transition rates for the flip of the spin 

W (S k -S fc ) = G (flfc)e~" S * B * (20) 
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where the function Gq(B) is an even positive function of B (see [25| for more details) 

G (B) = Go(-B) > (21) 
In the following, we will consider as two important examples the 'simple' dynamics 

Simple Dynamics : G s Q imple (B) = 1 (22) 

and the Glauber dynamics [23[ 

Glauber Dynamics : G9 lauber (B) = (23) 

v ' 2 cosh/35 v ' 

B. Mapping onto a quantum Hamiltonian 

As recalled in detail in [25[, the Master equation with the transition rates of Eq. [20] can be mapped via a similarity 
transformation onto the quantum Hamiltonian j26l - l30j 

H = E G ° ( E J (e-^(S^* J ^<) - of) (24) 

fc=l / 

for 2^ quantum spins described by Pauli matrices, with the following properties. The ground state energy exactly 



vanishes Eq = and represents the thermal equilibrium. The smallest non- vanishing energy E\(2 N ) determines the 
equilibrium time t eq (2 N ) (defined as the largest relaxation time of the Master Eq. [To]) 

C. Real Space Renormalization of the associated quantum Hamiltonian 

In our previous work [25], we have introduced a real-space renormalization procedure for quantum Hamiltonian of 
the form of Eq. [21]to determine E\ near zero-temperature. We now describe the application to the Dyson hierarchical 
model. The quantum Hamiltonian associated to the stochastic dynamics of the Dyson hierarchical Ising model can 
be rewritten as a sum of elementary operators (Eq I24[) 

2 jv-i 

H = + ( 26 ) 

B R ( -Per' Jj Qa * + EL\ \ B R 

h 2i -t + h 2l = G (J a* 2i + (e ^ 2 J - a^A + Go^o^ + -f] 

with B R is the static renormalized local field introduced Eq. [T5l 
Let us consider the first two spins : as explained in [25| . the sum 

hi+h 2 = Go(Jo(t| + —^-) ( e ~P <T *\ Jo<T % + ~2~ J _ r , j 1 _ : _ (; n (./ (l ,T- _1L) I , - ' ! - a\ \ (27) 



2 



K + y) e v '-of I 



can be renormalized onto the following operator describing the flip of the renormalized spin S R — (Si + 5 , 2)/2 (Eq 

□33 



h R _ r l 
"(1.2) = Ul 



?f (Bf ) (V^« s * - o%) (28) 
where the renormalized amplitude G R (B) can be computed from the function Bq via 



Gf{B) =e-V HB ™f°t §)G °_i J °-7} r \ (29) 



e' 3 fGo(Jo + f) + e-^Go(Jo-f) 
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or equivalently using the inverse 

1 _ 
G?(B) ~ 2 

In conclusion, the real space renormalization of the quantum Hamiltonian corresponds for the Dyson model to the 
renormalization of the function G(B) starting from the initial condition Gq(B) that defines the initial dynamics (Eq 

ED. 



Gq( Jq 



2 I 



Go (Jo 



B_ • 

2 - 



(30) 



D. Equilibrium time t eq (2 N ) for a finite system of 2 N spins 



We apply iteratively the previous renormalization rule as follows : 

(0) The initial dynamics concerns the Dyson hierarchical model of 2 N spins with the couplings (Jo, .., Jjv-i) and 
the function Gq(B) that defines the transition rates of Eql20l 

(1) After 1 RG step, we have a Dyson hierarchical model of 2 N ~ 1 renormalized spins with the couplings ( Ji, .., Jjv-i) 
and the function Gf (£?) obtained by Eg [BUI 



1 



Gf{B) 



E 



e f^ ± Go(Jo + eif ) 



(31) 



(2) After 2 RG steps, we have a Dyson hierarchical model of 2 N 2 renormalized spins with the couplings ( J2, .., Jjv— 1) 
and the function G^iB) obtained by Eq[30l 



VE 



G 2(B) 2 ^ ±Gl(Jl+e2 |) 



(32) 



(k) After k RG steps, we have a Dyson hierarchical model of 2 renormalized spins with the couplings 
(Jfc, .., Jjv-i) and the function Gf(B) 



E 



Gfc-l(Jfc-l + £fe 2") 



(33) 



(AT) After N RG steps, only one spin remains, whose flipping is governed by the function Gjy(B) 



G*{B) 



E 



Gjv-i(Jjv-i + -J ] 



e N =± 



(34) 



However since this spin is alone, there is no local field B = 0, so we just have to compute the final number G^ nal — 
G*(B = 0) 

1 1 „pj N -i 



G final G*(B = 0) G N -l{J N - X ) 



E 



-/3e f 



I JjV-2 + e iV-l 2 



E E 

cn — i=± ejv— 2=± 



2 e 



G 



iV-3 



Jjv -3 + £jv-2- 



■/iV-2+ejv-l - 



2 jv- 



— EE - E 



Go [So] 



(35) 
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in terms of the variables B n that can be computed from the recurrence 

B n 

B n -1 = Jn-1 + € n—T- 



(36) 



with the initial condition 



B N = 



The hrst terms read 

Bn-i = Jn-i 

Bn-2 = Jn-2 + Jn-i ( N y 1 ) 
and more generally for < n < N — 1, one obtains 

N-l 

B n = 



n — Jn + Jm ( 2 ) 



In the following, we will need 



m=n+l k—n+1 



N-l m 



Bo(ei) -,^n-i) = Jo + E J,n II (t 



m—1 k—1 



and the sum present in the exponential of Eq. 



S(ei, .., £n —i) = E en ~Y =Ey 

n—l n—1 
N-l 

= E Jm 
m—1 



N-l 



N-l 



J n + E ^ n (| 

m— n+1 fe= n+1 
m m 

En(?) 

_ri=l fe=n 



(37) 



(38) 



(39) 



(40) 



(41) 



In conclusion, the equilibrium time t eq (2 ) that can be computed from the final amplitude G ™ a 25 1 



c\s~i final 
ZLr N 



reads 



t eq (2 N ) 



E E -E 

ejv_ i=± eiv-2— i ei=± 



e -^S(ei,..,ejv-i) 

G (Bo(ei, ..,ejv-i)) 



(42) 



(43) 



where i3o( £ ij .., ejv-i) an d S(ei, .., ejv-i) are given in Eqs 1401 and 1411 To determine the leading behavior near zero 
temperature, we have to distinguish whether the initial function Gq(B) does not depend on /3 ( as in the simple 
dynamics of Eq |2"2"|) or depends on (3 ( as in the Glauber dynamics of Eq[23]). These two cases are studied respectively 
in the two following sections. 



IV. EQUILIBRIUM TIME FOR THE 'SIMPLE' DYNAMICS 
A. Real Space renormalization Solution 



For the initial condition of Eq. [22] that defines the 'simple' dynamics, the result of Eq|43]for the equilibrium time 
becomes 



^simple ^N ^ 



2^ 



E E -E e " 

ejv-l— i ejV-2=± ei=rt 



/3£(ei,..,ejv_i) 



(44) 
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where 



JV-l 



m— 1 



En(f) 



_n=l fc=n 



2 2 L 2 



2 

EAT-2 



1 + l 



£jV-3 , eAT-3 e AT-4 



2 2 

£iV-2£jV-3 

22 



ei£2 " 
2 2 . 

ejy-3 £ jV-4---ei " 

2^-3 
ew-2£Ar-3---ei ' 



+ 



)N-2 



(45) 



B. Leading term near zero temperature 

At low temperature, the sum of Eq|44]will be dominated by the minimal value of £(ei, .., ejv-i) of Eq. [41] Since 
£at_i appears only in the last line, and since the term between the brackets cannot change sign with respect to the 
first term 1, we are led to the choice 



e;v-i = -1 



(46) 



and we have now to minimize 

S(ei, .., ejv_2, ejv-i = — 1) = J\ + J2~^ 



e^v-2 
2 

Jn-i 



Jn-i 



1 + t 



T 63 



2 2 2 



£AT-3 . £Af-3£JV-4 



eAf-3eAr-4---£i 



2 iV-3 



(47) 



Since from Eq[2] 



Jat- 



>N-2 



- = J 2( 1 -^ N - 2 \l-2-°)>0 



(48) 



we are led to the choice 



e/v-2 



(49) 



and so on, so that the minimum is achieved when all e$ = — 1. The corresponding minimum reads using Eq. [5] 



S(ei = -l,e 2 = -1,..,£at_i 
iv-i 

3 



1 



N-l 



-1) = ^ J„ 



Ef-5 



.71=1 



m— n+1 



(50) 



so that finally 



^.simple ^2"^) 



"eg V" y — 

We may now use Eq[2]to compute explicitly 



e l 



iV-l 



Jo + 



'4-H) ? 



(51) 



ln[2 Jv t^ mp/e (2 Ar )] =/3 



N-l 



'4-nr 



= 13 Jo 



4 1 - 2 (1 - a ^ N -^ 1 1 - f-2- ff )( A '- 1 ) 

1 + -2 1 " CT + -2- a 1 '- 

3 l-2!- CT 3 1 + 2" CT 



(52) 
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This expression is valid for finite N and we may check the first terms 

ln[2^ e (2)] =f3J 

In [2 2 tf q mple (2 2 )] = [3 J [l + 3.2- CT ] = 

To obtain the leading behavior for large N, we should now specify the value of a 



3 

Jo + g J i 



(53) 



Case a < 1 : Power-law barrier 



For a < 1, Eq. 1521 becomes in terms of the length Ln = 2 N 



a<l : ]n[2 N t?™ ple (L N = 2 N )] = (3 J 



4 L 



3 2 1 



(l-a) 
N 



4 

1 + - 



f f 



3 2° 



1 3 2 ff + 1 



+ Q{L N °) (54) 



i.e. the energy barrier near zero temperature scales as the following power-law of the length Ln = 2 



N 



In [2 N t?™*> le (L N = 2 N )] ^ 4 J 



3(2 X - CT - f) 



-(l-a) 
J N 



(55) 



D. Case a — 1 : Logarithmic barrier 

For cr = 1, Eq. [52] becomes in terms of the length Ln = 2 N 



a = l : ln[2 iv ^ m f ie (i JV =2 JV )] = /3 J 

= i9Jb 



3 9 



0(L 



N i 



31n2 



ln(L 



N) 



+ 0{L- N l ) 



i.e. the energy barrier near zero temperature grows logarithmically with the length Ljv = 2 N 

ln [2 N t^ le {L N = 2 N )] 4 J 



31n2 



ln(J 



N 



(56) 



(57) 



E. Generalization to other dynamics where Go{B) does not depend explicitly on /3 

It is clear from Eq. [43] that for all dynamics where Gq{B) does not depend explicitly on (3, the equilibrium time 
will be dominated again by S m i„ near zero temperature, with the result (using Eq I4TH) 



t eq (2 N ) 



-/3S mi „ 



Go(jB (e = -l,..,cjv-i = -l)) 



Go yJo + Em=l 



(58) 



i.e. the dynamical barriers are the same as above in Eqs 1551 and Eq. [57] (the only difference is in the prefactor of the 
exponential). 



V. EQUILIBRIUM TIME FOR THE GLAUBER DYNAMICS 
A. Real Space Renormalization Solution 

For the initial condition of Eq. [231 describing the Glauber dynamics, Eq 1431 becomes 



Glauber /^N 
eq 



(2 iv ) 



1 e^E^o 1 ^ 



2G 



final 
N 



2 N 



E E -E e " 



■/3E(ei,..,e N -i) 



e /3Bo(eL,..,ejv- 1 ) + g-/3Bo( 



ei,..,ejv_ij 



59) 
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B. Leading term near zero temperature 



So at low temperature, we have to minimize over the = ±1 and over e = ±1 the combination (Eqs|45]and Eq I40[) 



N-l 



S(ei, ..,ejv_i) +eB ( 



En 

,n=l k=n 



+ e 



JV— 1 m 

Jo + E Jm n 

m=l fc=l 



J 0e + Ji|(l + e ) + J 2 | [l + |(l + e)]+J 3 |[l + | + ^(l + e) 



-Jn-2 

-Jn—i 



2 

e JV-l 



1 H 1 — h •■■ H —*f 5 (1 + e) 



1 



2 2 2 

£jV-2 £jV-2 £ Af-3 

2 2 2 



2 iV-3 

e Af-2€Af-3--- £ l 

2 N-2 



+ 



(60) 



The new factors containing e are not able to change the signs of the expressions between brackets, so we find again 
by recurrence (Eqs 1461 and [4"9"|) that this function is minimum by choosing all ej = — 1 and finally e = — 1. The 
corresponding minimum reads 



[E(ei, .., ejv-i) + eS (ei, .., ejv-i) 



E(ei = -1, .., e N -i = -1) - Bo(ei = -1, », ejv-i = -1) 



JV-l 



2J Jjy 



m—n+1 



N-l 



— Jo + J]] Jn 

m— 1 

1 



AT-l 



TO — 1 / -, \ fe 

2 



fe=i 



1 ■ ' \ 



m — 1 



We thus obtain the following leading behavior at low temperature (Eq l59|) 



±Glauber fnN 



(2 iv ) 



~ _L P /9[E^o Jm-[£(ei,..,ej V -i)+eB (e 1 ,..,e JV _i)] mi 1 

~ 2 W 

= J_ i^:;^+j.+^:,::/4i-(-i) m 1 

2 jv 



We may now use Eq[5]to compute explicitly 



ln[2 w ^ ouber (2 iV )] =(3 

= PJ 
= PJ 



2J + J " 2(1_ 



rn—1 
N-l 



4 






3 




2 



AT-l 



2 + 1 ^ 2 (i-<x)m - - 2 {1 - a)m (-- 



m — 1 



2 + _ 2 !- CT 2 1 ~ a 

3 1-2 1 "' 7 3 1+2-- 7 



This expression is valid for finite N and we may check the first terms 

In [2t^ auber (2)] = 2/3 J 
In [2 2 tf q auher (2 2 )] = J [2 + 2 1 -*] = [2 J + Ji] 

To obtain the leading behavior for large AT, we should now specify the value of a. 



(61) 



(62) 



(63) 



(64) 
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C. Case a < 1 : Power-law barrier 
For a < 1, Eg [63l becomes in terms of the system size Ln = 2 N 
<7<1 : ln[2 N tf q lauber (2 N )] ~ J 



32 
4 L 



32 



1-<T 



- 1 

(l-<r) 
N 



3 2 <T - 1 - 1 3 2 CT + 1 



l-CT 



„ 4 1 
2+ -- 



32 



cr-l 



1 3 2 CT + 1 



i.e. the energy barrier near zero temperature scales as the same power-law of the length Ln = 2 



N 



In [2 N t^ auber (L N = 2 N )] ^ 4J 

P ~ 3(2!- ff - 1) 



L 



(i- CT ) 

N 



(65) 



(66) 



and with the same prefactor as in Eq. 1551 even if the finite corrections of Eq [55] arc different the finite corrections of 
Eq [SJ] (more explanations on these finite differences are given in Appendix [XJ . 



D. Case a — 1 : Logarithmic barrier 



For cr = l, Eq |6"31 becomes in terms of the system size Ln = 2 N 



a = l : ln[2 N t% auber (2 N )] ~ /3J 

r a 

= 13 J 



2+-(N-l) 



2 




'4 


4" 


9_ 


= 13 Jo 


3 N ' 


" 9 



4 , r 4 

In iv at H — 

31n2 N 9 



i.e. the energy barrier near zero temperature grows with the same logarithmic barrier 

ln [2 N t% a » ber (L N = 2 N )] 4J 



31n2 



ln(L N ) 



(67) 



(68) 



with the same prefactor as in Eq [57] even if the finite corrections of Eq. [S5] and of Eq. [ST] are different (see more 
details in Appendix [A"f . 



E. Generalization to other dynamics where Gq(B) oc e 



-P\B\ 



From the analysis presented above, it is clear that dynamical barriers near zero temperature will remain the same 
for all dynamics where the amplitude Gq(B) displays the same exponential decay as the Glauber amplitude of Eq. [22] 



G (B) 



-P\B\ 



(69) 



As an example of other dynamics with the same exponential decay, we may cite the Metropolis dynamics corresponding 
to (see Eq 120)1 



-0\B\ 



(70) 



VI. 



CONCLUSION 



In this paper, we have studied the stochastic dynamics of the Dyson hierarchical one-dimensional Ising model of 
parameter < a < 1 via the Real Space Renormalization introduced in our previous work [25[ . We have shown that 
this renormalization procedure amounts to renormalize a single function G(B) that defines the transition rates of 
the renormalized dynamics. We have solved explicitly the RG flow for two types of dynamics, namely the 'simple' 
dynamics (and other equivalent dynamics of section TlV E[) and the Glauber dynamics (and other equivalent dynamics 
of section IV Ep . We have obtained that the leading diverging dynamical barrier is the same power-law ]nt eq (L) ~ 

j3 ( gpF^^iy J L 1_<T for a < 1 and the same logarithmic term : \nt eq (L) ~ [/3 (jj^) — l] lni for a = 1, even if finite 

corrections are different, as explained in Appendix [A"l 
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Appendix A: Link between dynamical barriers and the highest energy cost of a domain-wall 

1. Case of the Glauber dynamics 

For the Glauber dynamics, the dynamical barrier for a system of 2 N spins corresponds to the maximal energy cost 
of one domain wall inside the system 

-\n[2 N t^ r (2 N )] = max (u { */ N ~ k) - U? N S ) (Al) 

where U^n' 2 ^ represents the energy of the configuration where the first k spins are (—1), whereas all others spins 
are (+1). 

The property of Eq. lAll can be easily checked for small systems, for instance for the case N = 1 containing 2 1 = 2 
and the case N = 2 containing 2 2 = 4 spins fEq [64]) 

1 • r_ .Glauber/ yi _ j _ r ,(l,l) TT GS 



In [2t^ auber (2)] = 2J = U^- i} - U c 2 



1 



In [2H^ r {2 2 )} = 2 J + J x = U^ 3) - U2 S = - Uf s (A2) 

P 

For an arbitrary number N of generations containing 2 N spins, the correspondence of Eq. lAll is less straightforward 
because the energy cost of the configuration where the first k spins are (—1), whereas all others spins are (+1) reads 

JV-l 

jjiK^-k) _ v gs = 2 j oCo (fc) + Jj [2cj(k) + (1 - 2 Cj (k)) 2 1 --'S 3 _ 1 (fc)] (A3) 

i=i 

in terms of the coefficients Ci(k) G {0, 1} of the base-two decomposition 

JV-l 

k = ^( fc )2 l = c o(fc) + ci(*)2 + c 2 (fc)2 2 + ... (A4) 

i=0 

and of the corresponding partial sums for i = 0, .., N — 1 

i 

Ei(k) = j2cj(k)2> (A5) 
j=o 

Nevertheless, one can check that the energy barrier obtained in Eq. [62] corresponds to the energy cost of a domain 
wall located at fc» 

1 N—l . 1 / 1 \ m ~ 1 

iln(2^™-(2^)) =,2J +£j m =(c/ 2 ( ^' 2N - M -f/ 2 G « S ) (A6) 

^ m— 1 L _ 

where fc* is characterized by the following coefficients in the base- two decomposition IA4I 

C2p(fc*) = 1 

c 2 p+i(M = (A7) 

so that 

K = 1 + 2 2 + 2 4 + .. (A8) 



2. Case of the simple dynamics 

For the 'simple' dynamics, the correspondence of Eq. IA1I between the dynamical barrier and the maximal energy 
cost of a single domain wall does not hold, as can be seen already for the case N = 2 corresponding to 2 1 — 2 spins 
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and for the case N = 2 corresponding to 2 2 = 4 spins, since the dynamical barriers of Eq[ 



iln[2^ e (2)] =J 



/3 

iln[2 2 i^ e (2 2 )] =J + ^Ji 



(A9) 



are clearly different from the maximal energy cost of Eq IA2I This can be understood as follows on the case N = 1 
with two spins. The transitions rates associated to the simple dynamics (Eqs l20l and l22|) read 



]Y 3imple (++ -> H ) = W slmple (++ ^ h) = 

tf/ sim P ie (H > ++) =w simple {A > ) = 



,-PJo 



(A10) 



whereas the Glauber transition rates (Eg I23p are given by 



W 



Glauber 



(++ _> +-) = w Glauber ( 



-PJo 



e +f3Jo + e-/ 3 ^ 



TF 



Glauber 



(+- -> ++) = v^ G/auhe ''(+- -> — ) = 



3 +/3Jo 



1 _|_ g-2/3Jo 
1 



e +/3Jo 4- e ~l3Jo 1 _|_ e -2^J 



(All) 



For 2 spins, the equilibrium time is determined by the rate W(-\ — h — > H — ) to create a domain- wall when starting 
from one ground state (the time to eliminate the domain- wall is then negligible), and these two rates are respectively 
of order e~^ J ° for the simple dynamics and of order e~ 2 @ J ° for the Glauber dynamics, i.e. the dynamical barriers 
differ by a factor 2. One could argue that the Glauber dynamics (or other equivalent dynamics of section IV E|) is 
more 'physical', in the sense that all transitions rates remain bounded near zero-temperature, whereas in the 'simple' 
dynamics (or other equivalent dynamics of section HV E[) . transition rates corresponding to a decrease of the energy 
diverge near zero temperature. Nevertheless, we should stress that the difference between the dynamical barriers of 
the two dynamics remains of order O(l), whereas the leading terms found in the text for the case a < 1 (Eqs 1551 and 
[66]) and for the case a — 1 fEqs[57land[68f are the same. 
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